A Hybrid Interpolation Method for Fractional PDEs and Its Applications to Fractional Diffusion and Buckmaster Equations

This study presents a novel numerical method to solve PDEs with the fractional Caputo operator. In this method, we apply the Newton interpolation numerical scheme in Laplace space, and then, the solution is returned to real space through the inverse Laplace transform. The Newton polynomial provides good results as compared to the Lagrangian polynomial, which is used to construct the Adams–Bashforth method. This procedure is used to solve fractional Buckmaster and diffusion equations. Finally, a few numerical simulations are presented, ensuring that this strategy is highly stable and quickly converges to an exact solution.